General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 










r:'*' 








A:. 




I , S ecurity ,-li,ssilication of tilto, holly k,( nUMwet imd indcxiiuj onno 

y 1 ORIGINATING ACTIVITY fCor/»tJKi(o mithor) 


DOCUMENT CONTROL DATA - R & D 

( „lj„ir i,rl tmJ jiidexin,'. uunoinrion he cnterci/ when the overuti report is- clusKincd) 

|2n. REPORT SGCUrOTV C U A SGI P I C A TIOM 


The Ionosphere Research Laboratory 


2 b» GROUP 


|3. REPORT TITLE 


Global Exospheric Temperatures and Densities Under Active Solar Conditions 


•j 4. DESCRIPTIVE NOTES (Typo ol report .■md.iticlttisivc tintes) 

I Scientific Report 

s 5, AUTHORISI (First name, middle initiolp Inst name) 

I Bernard J. Wydra 

je. REPORT DATE 

j October 3, 19T5 

S aa.\ CONTRACT OR GRANT NO. 


distribution statement 


Supporting Agencies 


70. TOTAL NO. OF PAGES 7b. NO. OF REFS 

96 

90. ORIGINATOR’S REPORT NUMBER(S) 


PSU-IRL-SCI-436 


9b. OTHER REPORT NO<s> (Any other numbers that may be assigned 
this report) 


11. SUPPLEMENTARY NOTES 


I 12. SPONSORING MILITARY ACTIVITY 


The National Aeronautics and Space 
Administration 
The Office of Naval Research 


1 13. ABSTRACT / „ ^ r. A • 1 I 

Temperatures measured by the OGO-6 satellite using the 6300 A airglow | 

snectrum are compared with temperatures derived from total densities and N^ 
dLsities It is shown that while the variation of the total densities with latitude and 
i magAtic activity agree well with values used for CIRA (1972). the temperature 
i behavior is very different. While the temperatures derived from the N^ 

I were in much bLter agreement there were several important differences which 
i radically affect the pressure gradients. The variation of temperature with magnetic 
”tity ^ndTcated a Lasonal and local time effect and also a latitude and delay 
l-=mp variation different from previous density derived temperatures. A new 
magJtrindeA isTroposed thL is better correlated with 

The temperature variations at high latitudes were examined for three levels of 

magnetic activity for both solstices and equinox conditions. A temperature 

maximum hi the pre- midnight sector and a minimum in the noon sector were noted 
and seasonal and geomagnetic time and latitude effects 

ture, density, pressure and boundary oxygen variations for the great storm of 
I March 8, 1970 are presented. 


FORM 

1 NOV 6 5 


1473 


(PAGE 1) 


S/N OrOl -807-6811 


NONE 

Security Classification 


PSU-IRL-SCI-436 , , / i 

Classification Numbers 1.4,2, 1.5.3, 1.6.1, 

1.9.4, 3.2.2 


Scientific Report 436 

Global Exospheric Temperatures and 
Densities Under Active Solar Conditions 

by 

Bernard J, Wydra 


October 3, 19V5 

The research reported in this document has been by The 

Office of Naval Research under Contract Grant No. N00014-67 A 
0385-0017 and the National Aeronautics and Space Administration 
under Contract Grant No. NGL 39-009-003. 


Submitted by: 



L./A. Carpenter, Assistant 
Professor of Electrical Engineering 


Approved 

Ionosphere Research Laboratory 
The Pennsylvania State University 
University Park, Pennsylvania 16802 



J. S. Nisbet, Director 
Ionosphere Research Laboratory 



ACKNOWLEDGEMENTS 


The author wishes to express his appreciation to Drs. J. S. 
Nisbet and L. A„ Carpenter for their advice and guidance during the 
course of this work and to Dr„ J. R, Mentzer for serving on his 
committee. Also he wishes to extend his gratitude to Drs. C. A. 
Reber and A. E. Hedin of the Goddard Space Flight Center for their 
assistance in using the OGO-6 density data, Dr, J, E. Blamont of 
C.N.R.S. in France, for providing the OGO-6 temperature data. Bob 
Divany, Beatrice Romaniec, and Beverly Beiswenger for their assis- 
tance in the computer programming, and to Susan, Kristin, Ann and 
Ursula Wydra for their support and tolerance. 

This work was supported by the Office of Naval Research 
under Grant No, N00014- 67- A- 0385- 001 7 and the National Aeronautics 
and Space Administration under Grant No. NGL 39-009-003, 


n 


TABLE OF CONTENTS 




ACKNOWLEDGEMENTS „ „ o . « . » . • » c . . 

LIST OF TABLES „ „o o o o . 

LIST OF FIGURES . , o , c » « » » , » o , . = * = 

ABSTRAGT » ° ■ 

CHAPTER I INTRODUCTION .o n o = ..»•• • 

1„ 1 General Statement of the Problem , « o » « • 

lo2 Development of Thermospheric Models . » . « 

1.3 OGO-6 Satellite Neutral Temperatures and 

Densities . . o o . » ° • • » • • 

1„4 Density Variations During Geomagnetic Storms . 
1,5 Temperature Variations During Geomagnetic 

Storms „ „ , o o o . o • 0 0 » ' • _ " • 

lo6 Geomagnetic Storm Effects on Thermospheric 

Wind Systems , „,o. o <..•••<>» • 

CHAPTER II METHODS OF DATA ACQUISITION , o . 

2, 1 OGO-6 Satellite Measurements . , « . . o . 


CHAPTER III DISCUSSION OF ANALYSIS = . 



3, 1 

Comparison of the Temperatures Measured 
from the 6300 A Line Width v/ith Model Tempera- 



tures Derived from Molecular Nitrogen and 
Total Mass Density o , 


3,2 

Variations of Temperatures as a Function of 
Magnetic Index Kp 


3o3 

Fourier Analysis of Quiet Day Latitude 
Variation at Solstice and Equinox , . , , » • 


3,4 

Anomalous Temperature Behavior in the 
Southern Hemisphere ,. ,,..«•••« • 


3,5 

Temperatures at High Geomagnetic Latitudes, , 


3,6 

Description of the Temperature N^ Density, 
O Density, Pressure and Boundary Oxygen 
Density Variation During the Great Storm of 



8 March 1970 ° ' 

'i 

3,7 

Time Constants of Thermospheric Effects Due 
to Magnetic Activity 




iii 


Page 

ii 

V 

vii 

X 

1 

1 

2 

5 

8 

13 

15 

19 

19 

21 

21 

27 

33 

34 
45 

68 

79 


Page 


CHAPTER IV CONCLUSIONS 85 

4„ 1 Comparison of the 6300A Temperature 
with Temperatures Derived from Total 
Density and from N« Density 85 

4, 2 Temperature Variations as a Function of 

Kp Magnetic Index 86 

4.3 Equinox and Solstice Global Variations 

for Quiet Magnetic and Solar Conditions .... 87 

4.4 Anomalous Southern Hemisphere Tempera- 
ture Longitude Variation ...... ... 87 

4. 5 Temperature Variation as a Function of 

Geomagnetic Time and Geomagnetic Latitude , . 88 

4.6 Temperature,. Density, Pressure and 

Boundary Oxygen Behavior Before, During 


and After a Large Magnetic Storm ..... 88 

4.7 Temperature Relaxation Time Constant .... 89 

REFERENCES 91 


iv 


Table 


LIST OF TABLES 


Page 

29 


1 AT/ AKp Coefficients for Summer 1970 Afternoon « » 

2 AT/ AKp Coefficients for March 1970 Morning « « » 

3 AT/ AKp Coefficients for Fall 69, Spring 70 and 

Summer 70 with AM and PM Local Times » » <> • o 

4 Coefficients of Fourier Series Fitted to Quiet Day- 
Orbital Temperature Variations at Equinoxes and 
Solstices oOOOOOOOOOOO"»»«»» 

5 Exospheric Temperatures and Accuracy of Averaged 
Data Tabulated by Geomagnetic Latitude and Geo- 
magnetic Time for Summer with Kp Range of 0 to lo3 . 

6 Exospheric Temperatures and Accuracy of Averaged 
Data Tabulated by Geomagnetic Latitude and Geo- 
magnetic Time for Summer with Kp range of 1<,3 to 

4 o 0 „ a o o o o . = o » 0 0 » - " ' ° ° 

7 Exospheric Temperatures and Accuracy of Averaged 
Data Tabulated by Geomagnetic Latitude and Geo- 
magnetic Time for Summer with Kp Range of 4» 0 to 

6o7o0<.0»0000000».0«0000 

8 Exospheric Temperatures and Accuracy of Averaged 
Data Tabulated by Geomagnetic Latitude and Geo- 
magnetic Time for Equinox with Kp Range of 0 to 1 = 3 . 

9 Exospheric Temperatures and Accuracy of Averaged 

Data Tabulated by Geomagnetic Latitude and Geo- 
magnetic Time for Equinox with Kp Range of 1»3 to 

4 „ 0 „ . » o O » O o o . » O O O O 0 » 0 » » 

10 Exospheric Temperatures and Accuracy of Averaged 

Data Tabulated by Geomagnetic Latitude and Geo.- 
magnetic Time for Equinox with Kp Range of 4» 0 to 

6 o 7 o . 0 o 0 o O o 0 0 = O . O O . 0 o o O 

11 Exospheric Temperatures and Accuracy of Averaged 

Data Tabulated by Geomagnetic Latitude and Geo- 
magnetic Time for Winter with Kp Range of 0 to 1 o 3 • 

12 Exospheric Temperatures and Accuracy of Averaged 
Data Tabulated by Geomagnetic Latitude and Geo- 
magnetic Time for Winter with Kp Range of 1.3 to 


V 


31 

33 

40 I' 

I 

50 

51 ! 

52 

56 

57 

58 [ 

62 1 

1^. 

U 

63 B 

il 

t:: 

i ■ 

t 

i; 


Table 


Page 


13 Exospheric Temperatures and Accuracy of Averaged 
Data Tabulated by Geomagnetic Latitude and Geo- 
magnetic Time for Winter with Kp Range of 4. 0 to 


64 


LIST OF FIGURES 


Figure 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


Comparison of 6300 A temperatures with and total 
density derived temperatures for an undis turbed period. 

Comparison of 6300 A temperatures with and total 
density derived temperatures for a moderately 
disturbed period , , , , , , , . , , , « . , . 


Comparison of 6300 A temperatures with and total 
density derived temperatures for a greatly disturbed 
pGriod 000000000,00000, 


AT/ AKp versus magnetic latitude derived from 6300 A 
temperatures, and total density derived tempera- 
tures, Jacchia (1^71) and Roemer (1972), for 1970 
Summer PM local time , , . . , . . » , , . 


AT/ AKp versus magnetic latitude derived from 
6300A temperatures, N 2 and total density derived 
temperatures, Jacchia (1971) and Roemer (1972), for 
1970 March AM local time , . , . . . . , . , 

AT/ AKp versus magnetic latitude calculated from 
6300 A temperatures for Fall 69, Spring 70 and 
Summer 70 with AM and PM local times . , . . . 

Orbital 6300A temperature distribution with fourier 
series fit for June 69 solstice , . , , , o . , , 

Orbital 6300 A temperature distribution with fourier 
series fit for September 69 equinox 

Orbital 6300 A temperature distribution with fourier 
series fit for December 69 solstice , , , « . , • 

Orbital 6300 A temp^erature distribution with fourier 
series fit for March 70 equinox , 

Orbital 6300 A temperature distribution with fourier 
series fit for June 70 solstice . . . . , , . , , 

6300 A temperature data versus longitude for magnetic 
latitude between -50° and -40° . . . , . . , , . 


6300 A temperature data versus longitude for magnetic 
latitude between -40° and -30*^ , , . , , , . , = 

6300 A temperature data versus longitude for magnetic 
latitude between -3 0° and -20° , , . , , , . , . 


Page 

23 

25 

26 

28 

, 30 

. 32 
, 35 

• 36 

o 37 
0 38 
. 39 
. 41 
. 42 
. 43 


vu 


Page 


Figure 

16 

16 

17 

18 

19 

20 
21 
22 

■1 

23 

24 

25 

26 

27 


6300 A temperature data versus longitude for magnetic 
latitude between -20° and -10° o 


OOOOOOffOO 


Exospheric temperature as a function of geomagnetic 
latitude and geomagnetic time for Kp range of 0 to 1.3 
in Summer . . » » . « » » » » » » ° 


O 0 0 


Exospheric temperature as a function of geomagnetic 
latitude and geomagnetic time for Kp range of 1.3 to 
4o 0 in Summer . . . = o » ° ° ° ° 

Exospheric temperature as a function of geomagnetic 
latitude and geomagnetic time for Kp range of 4.0 to 
6. 7 in Summer 


0 o 0 


OOOOOOO -8 0000 


Exospheric temperature as a function of geornagnetic 
latitude and^ 'geomagnetic time for Kp range of 0 to 
1„3 at equinox „ . . . » . » » ° ° ° ° ° ° ° 

Exospheric temperature as a function of geomagnetic 
latitude and geomagnetic time for Kp range of 1 <, 3 to 
4. 0 at equinox . . . .. . o . « » » = ° ° 

Exospheric temperature as a function of geomagnetic 
latitude and geomagnetic time for Kp range of 4. 0 to 

6.7 at equinox „ . . . . . . . . ° ° ° 

Exospheric temperature as a function of geomagnetic 
latitude and geomagnetic time for Kp range of 0 to 

1 „ 3 in Y/inter . . . . . . » « « » “ ° ° ° 

Exospheric temperature as a function of geomagnetic 
latitude and geomagnetic time for Kp range of 1.3 to 
4. 0 in Winter ... o.. .o. . = .■>«» » 

Exospheric temperature as a function of geomagnetic 
latitude and geomagnetic time for Kp range of 4.0 to 
6. 7 in Winter 


0 O 0 0 8 o 0 


0 0 8 0 


o O O 8 o o o 


Indicence of auroral forms during 1964 and 1965 
after Stringer and Belon (1966) 

6300 A. temperature versus magnetic latitude for 
undisturbed, moderately disturbed and greatly 
disturbed periods in March 1970 


OOOOO OOOO 


44 

47 

48 

49 

53 

54 

55 

59 

60 

61 

67 

69 


6300 A temperature (in hundreds of °K) contour map 
representing the 300 km variations witn magnetic 
latitude and universal time for five days in March 
1970 ... o o . 0 » 0 ° ° ° 


70 


vui 


.1 


J . - . . I 


Figure 

_ 3 

28 LogjQ (N 2 number density (m )) contour map 
representing the 400 km variations with magnetic 
latitude and universal time for five days in 

March 19^0 „ooooooooi)ooo«»« 

_ 3 

29 Log, - (0 number density (m )) contour map 
representing the 400 km variations with magnetic 
latitude and universal time for five days in 

March 1970 „„„o,,ooo<.oooo<.o 

30 AE and Kp indices for the period of March 8 to 12, 

1970 ooooooooooooooo""”" 

31 Variation of pressure at 400 km versus magnetic 
latitude for undisturbed, moderately disturbed and 
greatly disturbed period in March 1970 „ 

32 Variation of lower boundary atomic oxygen versus 
magnetic latitude for undisturbed, moderately 
disturbed and greatly disturbed periods in March 1970 

33 Atomic oxygen at 400 km versus magnetic latitude 

for undisturbed, moderately disturbed and greatly 
disturbed periods, and 3 hours after the greatly 
disturbed period in March 1970, « <• , » » ° • « 

34 Standard deviation of least squares fit of 6300 A 
temperature data versus relaxation time constant 
T for Fall 69, Spring 70 and Summer 70 for AM 

and PM local times o» , = ° ° 

35 AT/ AW versus magnetic latitude for Summer AM 

and PM local times „ 0 , »o » 

36 Exospheric temperature versus magnetic latitude 

for W = 0 and 10,7 cm solar flux =150 for 
Summer AM and PM local times , , , , 0 , , . 


Page 

71 

72 

73 

75 

77 

78 

82 

83 

84 


IX 


j o ■ ; I i 


ABSTRACT 

Temperatures measured by the OGO-6 satellite using the 
6300 A airglow spectrum are compared with temperatures derived 
from total densities and N 2 densities. It is shown that while the 
variation of the total densities with latitude and magnetic activity 
agree well with values used for CIRA (1972), the temperature behavior 
is very different. While the temperatures derived from the N 2 
density were in much better agreement there were several important 
differences which radically affect the pressure gradients. The varia- 
tion of temperature with magnetic activity indicated a seasonal and 
local time effect and also a latitude and delay time variation different 
from previous density derived temperatures, A new magnetic index 
is proposed that is better correlated with the observed temperatures. 
The temperature variations at high latitudes were examined for three 
levels of magnetic activity for both solstices and equinox conditions, 

A temperature maximum in the pre-midnight sector and a minimum 
in the noon sector were noted and seasonal and geomagnetic time and 
latitude effects discussed. Neutral temperature, density, pressure 
and boundary oxygen variations for the great storm of March 8, 1970 
are presented. 
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CHAPTER I 


INTRODUCTION 

1 . 1 General Statement of the Problem 

Satellite measivrements of total density date from 1957 when the 
orbital decay of Sputnik I was used to calculate atmospheric densities. 
The total density varies greatly with altitude decreasing by a factor of 
10 in 100 km so that some method of unifying the data had to be found. 
The solution to the problem was provided by the work of Nicolet (1961) 
who developed a family of atmospheric models based on the tempera- 
ture profile assumed by the atmosphere when allowed to cool from a 
high temperature. By using constant boundary conditions these models 
enabled a unique relation to be established between the density at any 
altitude and the exospheric temperature for the appropriate model. 

Stein and Walker (1965) showed that depending on the boundary- 
conditions for densities, temperatures, and temperature gradients a 
wide variation in models could be used to fit the experimental data. 

The models of Jacchia (1964, 1965, 1970, 1971) used the total 
density measured by satellite drag to derive the exospheric tempera- 
ture, Hedin, et al. , (1972) used the molecular nitrogen density 
measured by a quadrupole mass spectrometer on the OGO-6 satellite 
to construct a model of the thermosphere, and to derive the exospheric 
temperature. Blamont and Luton (1972) used a Pabry-Perot 6300 A 
interferometer on the OGO-6 satellite to measure the thermospheric 
temperature at 250 to 350 km. 

Now, using both the neutral temperature and density data from 
the OGO-6 satellite it is possible to compare the different methods of 
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obtaining the exospheric temperatures with thermospheric tempera- 
tures measured at approximately the same time and place,, The 
combined data sets also makes it possible to compare the variations 
of temperature as a function of magnetic activity for each of the 
different methods of deriving the exospheric temperature, to relate 
them to previous measurements, and to explain the differences 
between themo 

1 . 2 Development of Thermospheric Models 

Nicolet (1961) examined satellite drag results and initiated a 
model of the thermosphere based on the theoretical calculations in 
which the atmosphere was allowed to cool. When this was done it was 
found that consistent temperature and density profiles were obtained 
which were almost independent of the initial conditions assumed,, 

Jacchia (1965, 1970, 1971) used satellite drag total density to 
construct a static diffusion model of the atmosphere above 110 km. 
The CIRA model incorporates the Jacchia (1971) model and also 
presents a review by Roemer of the new data which agrees and 
disagrees with the model, Jacchia (CIRA 1972) proposed that the 
model will show poorer results for comparatively short term varia- 
tions - such as diurnal variation or geomagnetic storms - than long 
term variations , As discussed by Bauer (1967), Nisbet (1967), 

Swartz and Nisbet (1971) and Bauer, Waldtuefal and Alcayde (1970), 
incoherent scatter radar neutral temperature measurements indicate 
a phase difference of about two hours between the Jacchia (1971) 


model and measured data. 
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An einpirical model based on OGO-6 mass spectrometer data 
has been constructed by Hedin, et aL , (1974) for geomagnetically 
quiet conditions,, This model unlike CIRA (1972) includes the morning 
maximum for helium, the high latitude of the diurnal temperature for 
solstice conditions while the oxygen maximum remains at low latitudes 
and the helium maximizes in the winter hemisphere, a winter 
maximum in the O to N£ density ratio at 120 km with a 400 K summer 
to winter temperature difference, and the 1600 hour maximum for 
with an extended bulge in latitudeo However, Hedin, et alo , (1974) 
suggested cautious use of the model results for conditions which are 
beyond the data base used in generating the model. 

Chandra and Herman ( 1969 ) solved the equations of continuity 
and heat conduction simultaneously for electrons, ions and neutral 
constituents. The model predicts a decrease of atomic oxygen in the 
lower thermosphere which results in a decrease of the peak F2 
electron density and an increase in the neutral temperature and 
density. At higher altitudes the O density increases which causes the 
electron density to increase. This means that the slab thickness of 
the F-region increases, though the peak electron density and total 
electron content decrease as has been reported by Evans (1970), The 
and 0„ at higher altitudes should also increase as a result of the 

Cd Ct 

higher neutral temperature. 

Perturbation theory formed the basis for a three-dimensional 
model by Volland and Mayr (1971) which was used to analyze a nights 
time source of heat input in the auroral zone. The results indicated 
that the effect of the disturbance depended on storm time, local time, 
and individual longitude of each event. The calculated delay, time 


varied from 2 to 8 hours from high to low latitudes,, 

The thermospheric wind circulation was investigated by Mayr 
and Volland (1972) using Joule heating as the energy source. The 
results indicated that the molecular nitrogen expansion in the auroral 
zone generates meridional winds towards the equator. The atomic 
oxygen density variation at 300 km with increasing magnetic activity 
is negligible while at 400 km it increased slightly. The molecular 
nitrogen density at high latitudes is ten times the concentration at low 
latitudes for Ap = 200, The model was used to reproduce and explain 
OGO-6 mass spectrometer data, Mayr and Volland (1973) used 
another three-dimensional two constituent model to describe the 
effects of geomagnetic activity. The model exhibits steep tempera- 
ture enhancement at high latitudes, the depletion of helium and atomic 
oxygen, the negligible latitude variation of atomic oxygen at the same 
time that molecular nitrogen exhibits a gradual but significant increase 
towards the auroral zone, and the change of delay time from 2 to 8 
hours from high to low latitudes. The results compared favorably 
with OGO-6 mass spectrometer data, Mayr and Volland (1974) 
proposed that adiabatic heating by thermospheric wind circulation 
would increase the neutral temperature at low latitudes. Calculations 
showed that the winds generated in the auroral region would peak in 
less than two hours at midlatitudes after the peak of the heat input, 
Mayr and Volland (1974) concluded that OGO-6 observations can be 
explained by diffusive mass transport associated with the thermo- 
spheric circulation, 

Reber, Hedin and Chandra (1973) further analyzed the OGO-6 
data and concluded that dynamics play ii significant role in the upper 
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atmosphereo Specifically, the mass spectrometer data indicates that 
during geomagnetic disturbances circulation cells are formed which 
decrease the atomic oxygen at low altitudes where the lower boundary 
of the static diffusion models have been established. The phenomena 
observed by OGO-6 substantiates the need for dynamic processes to be 
included in the thermospheric models if the models are to reproduce 
the variations of density and neutral temperature during geomagnetic 
distrubanceSo 

1^3 OGQ-6 Satellite Neutral Temperatures and Densities^ 

A Fabry- Perot interferometer experiment to measure the 
neutral temperature at altitude around 250 to 320 km is described by 
Blamont and Luton (1972). The OGO-6 satellite was launched into an 
82° inclination polar orbit so that the local time observed by the 
satellite changed slowlyo The 6300 A airglow measurements enabled 
Blamont, Luton and Nisbet (1974) to determine the global temperature 
distribution for geomagnetically quiet conditions and each of the 
equinoxes and solstices. The conclusions of the quiet condition study 
are that the maximum temperature occurs near the summer pole at 
the solstices and that the effect of the day length at high latitudes is 
very pronounced, 

A quadrupole mass spectrometer on the satellite provided an in 
situ measurement of atomic oxygen, molecular nitrogen and helium 
densities, Hedin, et al, , (1974) used the undisturbed condition data to 
form a global model of the atmosphere at 450 km and 120 km. The 
model agreed well with both static diffusion models and incoherent 
scatter radar results, Reber and Hedin (1974) used the density 
data to propose a persistent heat source in the high and middle 
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latitudes during magnetically quiet periods. The heat source 
appeared to be co- rotating with the magnetic poles, exhibiting a 
maximum around 900 and 2100 hours universal time. The effect is 
present throughout the year and is stronger at the equinoxes than in 
the summer, and may have a maximum during the winter solstice. 

Three studies of thermospheric response to geomagnetic 
disturbances using OGO-6 satellite data were reported. Reber, 

Hedin and Chandra (1973) deduced exospheric temperatures from the 
molecular nitrogen densities to demonstrate that the temperature at 
middle latitudes increased more than 1000°K while the low latitude 
temperature increased by only 50° to 150°K during the great storm of 
8 March 1970. The authors also suggest that the energy is deposited 
in the high latitudes during the storm, and thereby raises the tempera- 
ture of the thermosphere and exosphere and generates a large scale 
circulation cell in the lower thermosphere, 

Blamont and Luton (1972) used the 6300 A airglow related 
neutral temperatures to describe the storm period between September 
25 and October 5, 1969. The results substantiate the dependence of 
the geomagnetic perturbation on latitude, but unexpected maximums 
of temperature occurred at the polar regions during quiet times. The 
observations shov/ed that the principal effect of the geomagnetic 
activity occurs in polar regions and appears to depend on local time. 
The zone where this strong effect occurs is at lower latitudes at night. 
Very high temperatures are reached in this zone compared to the day 
and this may result in strong temperature gradients in the meridional 
direction which induce neutral winds with velocity above 500 m/ sec. 
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The features of the storm temperatures demonstrate that the high 
latitudes are quickly affected by geomagnetic activity with a delay of 
less than 3 hours and the equatorial region has a more gradual tem- 
perature variation with a smaller magnitude,, 

Taeusch, Carignan and Reber (1971a, 1971b) also investigated 
the period between September 27 and October 3, 1969. The observa- 
tions confirm many of the features reported by Blamont and Luton 
(1972). The density above 400 km reacted very quickly (less than an 
hour) to the geomagnetic activity at high latitudes, but was relatively 
unaffected in the equatorial region. The O/N^ ratio variations 
suggested dynamic processes which result in a circulation that is up- 
ward at the poles, then towards the equator where the process 

subsided. Taeusch, etal., (1971a, 1 97 lb) concluded that the major 
effect of the geomagnetic storms during the observation period was in 
the high latitudes (above 50° magnetic latitude) and that the data is 
better correlated with magnetic latitude than geographic latitude. 

Anderson (1973) used a microphone total density gauge on the 
OGO-6 satellite and also investigated the storm period of September 
25 to October 3, 1969 for magnetic latitudes between the magnetic 
equator and 40°N. A correlation study of the densities and magnetic 
indices indicated that the low latitude density response lagged the Dst 
index very little, while the density response lagged the Ap index by 
about 3 hours. The magnitude of the increase in total density was 
different from that suggested by Taeusch, et al. , (1971a, 1971b). 
Anderson reported that the total density was over twice as large for 
Kp = 8 as for Kp =0. 
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Many investigations have used various indices to try to quantify 
the variations of temperature and density with respect to the intensity 
of geomagnetic storms. An excellent review of the derivations of the 
Kp, Ap, AE, and Dst indices has been written by Rostoker (1972) in 
which he also discusses the strengths and weaknesses of each index. 

The main general weakness of all the indices is the non-uniform 
distribution of stations used to obtain the global indices. Small 
spatial variations with small time durations may be entirely missed by 
the measuring stations and therefore invalid conclusions drawn from 
the reported. indices , 

1 ^ 4 Density Variations During Geom agnetic Storms 

Jacchia (1959a, 1959b) analyzed the speed and attitude changes 
of spherical satellites and deduced the total density necessary to 
produce the observed effects. The analysis of satellite drag was a 
major breakthrough in the development of the global picture of density 
variations. Diurnal, seasonal, latitudinal and longitudinal changes 
were able to be uniformly described on a global scale, Jacchia (1965, 
1970, 1971) formulated a model based on the satellite drag density 
measurements which he has periodically reviewed and revised to 

Include new data (Jacchia, 1974), 

Much of the new density data analyzed has been used to attempt 

a description of the effects of geomagnetic activity on the global 
density structure. Lew ( 1969 ) used tow altitude satellite drag data to 
demonstrate that the time delay between the maximum magnetic 
activity and the maximum variation of the densities below 200 km was 
between 2,4 and 11,52 hours, depending on the intensity of the storm. 
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Lew (1969) used the planetary Kp as the magnetic activity indicator, 
and analyzed the data from two different satellites „ However, he was 
not able to determine the latitude dependence of the delay time„ 

Roemer (1971a, 1971b) reviewed the observational results up 
to 1971 and ‘presented Satellite drag measurements from six satellites 
to reach several conclusions„ The temperature inferred from the 
density varied from the quiet condition as a function of latitude given 
by 

AT {(py = [ (21„4 ± 8„3) sin j 0 | + 17» 9] Kp + 0„03 exp (Kp) 

( 1 ) 

where (j) is the latitude and Kp is the 0.4 day mean of the planetary 
index Kp. The mean time lag of the densities below 200 km is about 
12 hours whereas the time lag above 250 km is approximately 5. 5 ± 
0.3 hours after a geomagnetic disturbance. There was a dependence 
of time lag on the length of the storm which is given by 

At = (0. 13 ± 0,01) A T + 0.01 (2) 

where T is the storm length and time lag is given in days. The time 
lag above 250 km was not found to depend on altitude, local time or 
latitude. However, Reomer reported that the specific heating depend- 
ence is 30% higher at night which led him to believe that the heat 
source was either stronger or the thermosphere more sensitive at 
night. Roemer and Lay (1972) used the same six satellites, orbiting 
from 1961 to 1966, to reiterate the time delay independence of 
latitude, and that the magnitude of the variations did depend on 
latitude and longitude. 
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A technique similar to satellite drag analysis has been used by 
May and Miller (1971) in which the spin rate of satellites equipped 
with solar cell paddles was observed„ The Ariel 2 spin rate derived 
density showed a better correlation of low latitude variations with Dst 
than with Ap, but this may be due to the fact that Dst is measured at 
equatorial stations, whereas Ap is measured at high and mid- latitudes „ 
In any case, the densities exhibited a time lag of about 2 hours with 
the Dst and an increase of 1 gamma in Dst was equivalent to a 1 „ 25 
percent increase in densityo The data also displayed a slow recovery 
of the density after the storm, to values measured before the storm„ 

Various methods of measuring density have been used since the 
early studies done with satellite drag. Reviews and comparisons 
have benn written by Champion (1969> 1972), Von Zahn (1970) and 
Roemer (in the Cospar International Reference Atmosphere, 1972) » 

Newton (1970) used data from the Explorer 32 satellite density 
gauges to describe a strong positive density gradient above 40° 
latitude which indicates a permanent source of heating above that 
latitude. He also observed that the heating may be influenced by 
increased geomagnetic activity, Anderson (1973) used a microphone 
density gauge on the OGO-6 satellite to measure total density, and his 
results exhibited a lag of about 3 hours near the auroral zone. The 
data showed that the density for Kp = 8 was twice the density for Kp = 0 
at about 20 degrees latitude, 

A low g accelerometer calibration system (LOGACS) was used 
by DeVries (1972a, 1972b) to measure total densitites during a period 
of high geomagnetic activity in May 1967, The time lag of the density 
variations in the polar regions was 90 minutes or less and gradually 
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increased to 4. 5 to 6 hours at the equator. This led him to conclude 
that most of the energy associated with geomagnetic activity is depos- 
ited in the auroral zones, and then transported toward the equator by 
atmospheric waves and convective circulation. The density values 
before the disturbance were larger in the high latitudes, but after the 
storm, the values of density at about 150 km in the sunlit auroral 
region were only 70% of those before the storm. At the same time the 
densities at around 15° latitude were double the values before the 
storm, so that the 150 km densities were higher at low latitudes after 
the storm. However, the density variation at 170 km did not exhibit 
the depletion phenomena, but remained higher at the high latitudes. 

The same data was further analyzed by Allen and Cook (1974) 
who concluded that there were three heating zones, near the equator, 
in the particle precipitation zones, and over the polar caps. They 
also suggested that particle precipitation and not Joule heating was 
the major source of heating. It was noted that the equatorial variation 
which was observed corresponded very well with the results of May 
and Miller (1971), and the heating in this region was probably due 
to the decay of the ring current. 

An accelerometer was used on the OVl-15 (Spades) satellite 
provided data which led Forbes and Marcos (1973) to conclude that 
enhanced density variations at high latitudes are associated with 
auroral electrojet activity as indicated by the AE indices. They 
support the Joule heating mechanism as the source of heating and 
also note wavelike structures in the density profiles which may 
indicate gravity waves generated in the auroral zone, Forbes and 
Marcos also suggest that anyone attempting to model the high latitudes 
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should incorporate electric fields, Joule heating and dynamics into 
the model. Truttse (1969), using 6300A emission measurements to 
deduce densities also noted a linear relation between density varia- 
tion and the ten hour sum of the AE indices. 

A radio frequency monopole mass spectrometer on the Cosmos 
274 satellite enabled Romanovsky and Katyushina (1974) to demonr- 
strate a morning- evening asymetry. , The morning density changes 
were significantly larger than the evening variations which led them 
to conclude that the effect is produced by precipitation of 10 to 40 Kev 
electrons, which has been shown to have a maximum in the morning 
hours. Prolss and Von Zahn (1974a, 1974b) used a similar gas 
analyzer on the ESRO-4 satellite to observe that the densities of 
molecular nitrogen can increase by 150% and atomic oxygen can 
increase by 50% in the mid and low latitudes during geomagnetic 
storms. They also proposed that the neutral upper atmospheric 
storms and ionospheric storms are coupled and complimentary effects 
of geomagnetic storms, and that the data concerning these phenomena 
could be used to supplement each other. Philbrick (1974) used data 
obtained from mass spectrometers on OV3-6 and OVl-15 satellites 
to demonstrate that molecular nitrogen has a large variability in the 
auroral zone, presumably due to energy sources located there. He 
also notes that the measurements exhibit significantly larger varia- 
tions than presently allowed by static diffusion models and that these 
models may be inadequate to describe the atmosphere during severe 
geomagnetic storms. 
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1 , 5 Temperature Variations During Geomagnetic Storms 

A method of determining neutral temperatures from para- 
meters measured hy incoherent scatter radar was first developed by 
Nisbet (1967) using the atomic oxygen thermal energy balance equation 
and data taken at the Arecibo radar site. The results show a signi- 
ficant difference in the phase of the diurnal variation of the tempera- 
ture and the diurnal variation of the satellite drag densities. There- 
fore, there is a discrepancy between the incoherent scatter radar 
temperature measurements and the Jacchia (1965, 1970, 1971) model 
temperatures derived from satellite drag densities (Swartz and Nisbet, 
1971), 

Bauer, Waldteufel and Alcayde (1970) used (St. Santin-Nancay) 
radar measurements of electron density, ion and electron tempera- 
ture to determine oxygen concentration and exospheric temperature. 
The results exhibited a similar phase difference between density and 
temperature to that reported by Nisbet (1967), and also found that the 
temperature shape parameter displayed higher values in the day than 
at night. 

Analytic Green's functions enabled Thomas and Ching (1969) to 

\ 

solve the heat conduction equations and to calculate the response of 
the thermosphere to time dependent heat sources. Their results 
showed that the time lags of temperature and density were distinctly 
different and that there was no simple relationship between them. 

The authors stated that the possible reason for this was that the 
density lag depended on pressure and the entire temperature profile, 
whereas the temperature lag depended only on the local temperature. 
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Magnetic storm effects on thermospheric temperatures were 
studied by Waldteufel and Cogger (1971) using data from the incoherent 
scatter radar at Arecibo, They proposed that storm effects could not 
be described by a simple relation of the magnitude of the planetary Kp 
index, and that the characteristics as well as the time delay of the 
observed effects vary for different storms. 

Truttse and Yurchenko (1971) used a Fabry-Perot interfero- 
meter at the Zvenigorod Observatory in the Soviet Union to measure 
the width of the Doppler contour of the atomic oxygen red line (6300 A) 
emission. The method to calculate the temperature from the measure- 
ment was similar to that of Blamont and Luton (1972) and the results 
showed higher temperature variations than predicted by static 
diffusion models ?jased on satellite drag density measurements 
(Jacchia, 1965, 1970, 1971), and that the temperature increased to- 
wards the north which would be in agreement with the OGO-6 results 
(Blamont and Luton, 1972, and Taeusch, Carignan and Reber, 1971a, 
1971b). 

Rees (1971) used fluorescent emission measurements made 
during rocket flights in the morning and evening twilight to demon- 
strate that the heating associated with geomagnetic activity was 
enhanced by a factor of 2 compared to the global mean Kp variation 
derived temperatures using the formula of Roemer (1970). Rees (1971) 
also suggested that the time required for the neutral atmosphere to 
relax to the undisturbed condition was about 12 to 24 hours after the 
end of the magnetic activity. 

Watkins and Banks (1974) applied a method similar to Nisbet 
(1967) to determine neutral temperatures at the Chatanika incoherent 
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scatter radar site using electron density and temperature, ion 

temperature and velocity measurements. The preliminary results 

indicate that during one storm in the summer of 1972, the neutral 

temperature exhibited no substantial increase from the undisturbed 

day or the mildly disturbed day„ But the authors pointed out that the 

2 

Joule heating or the most disturbed day was less than 5 ergs/cm sec. 
Solar EUV attenuation was measured by Knight, Uribe and 
Woodgate (1973) at sunrise and sunset from the OSO-6 satellite to 
determine the neutral temperature. Their results show that the 
geomagnetic storm effects are about twice as large at dawn as at 
dusk, using the Kp index as an indicator of geomagnetic activity. The 
operational time period of OSO-6 overlapped the operational period 
of OGO-6 although the OSO-6 satellite was in an equatorial orbit. The 
local time effect is similar to that which Blamont and Luton (1972) 
reported, Knight, Uribe and Woodgate also suggested that the Kp 
index was not a good indicator of short term geomagnetic activity, 

1 . 6 Geomagnetic Storm Effects on Thermospheric Wind Systems 
Johnson and Gottlieb (1973) described a wind system by using 
the latitudinal variation of oxygen and calculating the circulation 
system that would be necessary to correspond to the variations. It is 
a system constructed for low magnetic activity. The atomic oxygen 
is produced at the subsolar point by photodissociation of molecular 
oxygen, and then transported by the wind system toward the poles at 
equinox and the winter pole at solstice, Blum and Harris (1972) 
calculated a quiet condition wind system using the Navier-Stokes 
equation combined with an atmospheric model of Jacchia (1965) and an 


I 


1 . 


J 


16 

ionospheric model of Nisbet (1971), the Penn State Mk I model. The 
bium and Harris model also predicts that the winds travel from sub- 
solar point to antisolar point. Kohl and King (1967) used Jacchia's 
(1965) model and an assumed function of ion density to produce 
similar results, 

Burge, et al,, (1973) used the temperature latitudinal depend- 
ence on Kp derived by Roemer (1971a, 1971b) to determine a wind 
system at storm times. The results showed radical alterations of the 
quiet wind system. The preferential heating of the high latitudes 
during geomagnetic disturbance is sufficient to reverse the normal 
poleward daytime wind during equinox and winter conditions, Johnson 
(1974) combined available data from various sources concerning 
density variations and concluded that although during undisturbed 
periods the winds are towards the poles, when geomagnetic storms 
occur, the heat input at the auroral region may be equal to or greater 
than the solar heating at low latitudes so that there will be large 
scale winds away from the auroral region. 

Cole (1971) calculated that Joule heating could result in winds 
of around 500 m/ sec equatorwards from the auroral region, Chang, 
et al, , (1974) used magnetohydrodynamic equations to calculate winds 
resulting from Joule heating. Winds of the order of 1000 m/sec were 
obtained, but the model neglected viscosity which would significantly 
reduce the winds (Zimmerman and Rosenberg, 1972), 

A theoretical study of the thermosphere using Joule heating 
and particle precipitation as the energy sources was completed by 
Hays, etal, , (1973), They calculated vertical winds in the auroral 
zone of a few m/ sec which would in turn drive horizontal winds of a 


17 


larger scale in a cellular motion. The authors also suggested that the 
downward movement at lower latitudes would cause compress ional 
heating and a decrease in the mean molecular weight. 

From the O/N^ ratio variations measured by the mass 
spectrometer on OGO- 6, Taeusch, Carignan and Reber ( 197 la) 
indicated the observance of a reversal in the quiet time winter hemi- 
sphere wind pattern during a geomagnetic storm in 1969. DeVries 
(1972b) indicated that the LOGACS data also exhibited large density 
gradients which would correspond to considerable winds directed away 
from the auroral zone, contrary to the undisturbed circulation, 

A ground based Fabry-Perot interferometer measuring the 
Doppler shift of the 6300 A emission of oxygen was used by Hays and 
Roble (1971) to obtain wind measurement data during geomagnetic 
storms. The data taken in the auroral oval at 42° north latitude 
indicated wind speeds of 250 to 300 m/sec and 350 to 400 m/sec in two 
storm events during the night at an altitude of about 400 km, 

Vasseur (1969) used St, Santin-Nancay incoherent scatter 
radar measurements of ion drift velocity to infer quiet condition 
neutral winds during winter 1966 and summer 1967. The results were 
compared to theoretical models and showed some agreement and some 
discrepancies. Amayenc and Vasseur (1972) concluded that electric 
fields must be considered in the calculation of the neutral winds from 
the ion drift velocities. Results from February and July 1969 were 
used to show the better agreement with theoretical models, Salah and 
Holt (1974) presented similar results using Millstone Hill incoherent 
scatter radar measurements. 
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Evans (1973) inferred neutral winds for quiet and disturbed 
conditions in October 1968 using the incoherent scatter radar at Mill- 
stone Hill Observatory. The quiet condition neutral winds showed 
good agreement with theoretical models, but the magnetically disturbed 
conditions showed the neutral wind direction at midday to be toward 
the equator, which would be a reversal of the undisturbed wind pattern. 
Carpenter and Kirchhoff (1974) used Millstone Hill radar measure- 
ments to infer neutral winds for several periods, quiet and disturbed. 
One disturbed period in May 1972 exhibited a neutral wind daytime 
reversal from quiet condition patterns similar to Evans (1973). 
Carpenter and Kirchhoff (1974) included the effect of electric fields in 
the periods measured, 

Brekke, Doupnik and Banks (1973) used the incoherent scatter 
at Chatanika to measure ion drift velocities and derive neutral 
velocities. The results agreed fairly well with theoretical model 
patterns, but the authors suggested that there may be distortions to 
the calculated patterns due to the drag of the high speed ion flows in 
the higher latitudes, Brekke, Doupnik and Banks (1974) stated that the 
position of the radar with respect to the auroral oval during the 
measurement is very important in the neutral wind and electric field 
calculations. The results of measurements taken during a geomag- 
netically disturbed period, August 3-9, 1972, showed that the wind 
pattern was altered from quiet conditions as a result of local pressure 
gradients due to a heat source in the auroral oval. The equatorward 
motion at night was enhanced and the poleward motion in the afternoon 
was diminished. 
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CHAPTER II 

METHODS OF DATA ACQUISITION 

2, 1 GGO-6 Satellite Measurements 

The OGO-6 satellite was launched 5 June 1969 into an 82° 
inclination polar orbit with a perigee of about 398 km and an apogee of 
about 1100 km. The orbit was such that the local time at the equator 
decreased by eight minutes each universal time day. The change in 
longitude was -25° for two consecutive orbits and the time required 
for each orbit was about 100 minutes. 

The 6300 A oxygen emission line profile was determined by 
using a spherical Fabry-Perot interferometer aimed towards the 
horizon. B lament and Luton (1972) discuss the instrument, procedure 
of operation, raw data handling, and determination of the neutral 
temperatures of the experiment in detail. The accuracy of the 
measurement of the neutral temperatures determined by the profile of 
the 6300 A line was calculated to be ±65°K for an individual measure- 
ment, The spatial resolution of the measurement was approximately 
1° longitude, 6° latitude, and 30 km in altitude at the equator. Only 
the reduced data for the daylight period has been used because of the 
sensitivity of the measurement. The magnetically quiet period neu- 
tral temperatures have been examined and presented by Blamont, 
Luton and Nisbet (1974), 

The quadrupole mass spectrometer used in situ to measure the 
concentrations of O, He and N^, is described by Carignan and Pinkus 
(1968) and Hedin, etal,, (1972, 1973, 1974). The accuracies of the 
measurements which these authors have determined vary with 
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constituent and altitude„ Systematic uncertainties in the measure- 
ments of O and increase with altitude and are usually less than 
10% to 15% near perigee c No data was used where the estimated 
uncertainty was greater than 25% for N 2 and O and 50% for He» 

There is also a laboratory calibration uncertainty of about 10% to 15 %« 
The spatial resolution of the measurement is less than 1/2° latitude 
and longitude and less than 5 km in altitude-. Magnetically qiiiet condi- 
tion data have been analyzed and presented by Hedin, et aU , (1972, 
1973 , 1974 ) and Reber and Hedin (1974)„ Also Taeusch, Carignan and 
Reber (1971a, 1971b) presented the densities observed during the 
magnetically disturbed period near the end of September and beginning 
of October of 1969, which was the same period examined by Blamont 
and Luton (1972) using the 6300 A derived neutial temperature data. 


CHAPTER III 


DISCUSSION OF ANALYSIS 

3. 1 Comparison of the Te mperatures Measured from the 6300 A 

Line Width with Model Temperatures Derived from Molecular 

Nitrogen and Total Mass Density 

The use of 6300 A emission spectrum to determine the thermo- 
spheric temperature is a relatively direct method of kinetic tempera- 
ture measurement. The basic assumption is that the radiating oxygen 
atoms are in equilibrium with the ambient neutral atmosphere. The 
altitude scanned by the interferometer varied from 250 km to 350 km 
with an average value of about 285 km. Temperatures in this altitude 
range are normally within 100 K of the exospheric temperature, A 
correction was applied to convert the temperatures to exospheric 
temperature using the Bates (1959) formula to facilitate comparison 
with other data 

T(Z) = T^ - (T^- T^) exp (-S(Z - Z^)) (3) 

where T(Z) is the meas'ared temperature at altitude Z and Z^ = 120 km, 

T = 355 K and S = 0.016 km" 
o 

To convert neutral densities to exospheric temperatures 
requires assumptions about the lower boundary conditions, the shape 
of the temperature profile and diffusive equilibrium of the atmospheric 
constituents. In this work it was intended to compare our values with 
the values of temperature derived from techniques adopted previously 
by others, and procedures similar to those that had been employed 
were adopted. , 
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The CIRA (1972) reference atmosphere and the model due to 
Jacchia (1971) were based on satellite orbital drag data. In the 
present work total density derived temperatures were taken to be the 
exospheric temperature of the appropriate CIRA (1972) model atmo- 
sphere that would give the measured total density at the altitude of the 

mea sur ement . 

The Goddard Space Flight Center OGO-6 model of the quiet 
thermosphere (Hedin, et ah, 1972) based its temperatures on mole- 
cular nitrogen densities. Molecular nitrogen is the major constituent 
in the region where eddy diffusion effects predominate and the effects 
of diffusive equilibrium assumptions are small. Boundary densities 
and the temperature profile still have to be assumed. The technique 
employed in this work was to match the observed N 2 density to the 
appropriate CIRA (1972) model. Temperatures derived in this 
manner are referred to as density derived temperatures." 

Figure 1 shows the magnetic latitude variations of the 6300 A 
temperature measurements and the total density and density 
derived temperatures during 11 March 1970 where the Kp magnetic 
index had an average value of 0.58 and a maximum value of 1.3. It 
is apparent that the total density and density derived temperatures 
are in excellent agreement at all latitudes, but that the 6300 A 
temperatures are systematically lower by 100 K to 150 K at lower 
latitudes. These observations are in agreement with Hedin, et al„ , 
(1972) who showed that the temperature estimated from their model 
was about 100°K larger than those observed by incoherent scatter 
techniques at Millstone Hill Observatory at 42. 6°N by Salah and 


Evans (1972). 
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The latitudinal temperature difference from the 6300 A airglow 
between equator and pole under these quiet conditions is observed to 
be about 300°K compared to the 150°K indicated by the temperatures 
derived from the densities „ 

Figure 2 shows the magnetic latitude variation of the 6300 A 
temperature measurements and the and total density derived 
temperatures for 6 and 7 March 1970 when the average value of the Kp 
magnetic index was 4.04 and the maximum value was 5. 7. Below 35° 
magnetic latitude the 63 00 A temperature is stilt about 150°K lower 
than the total density and N 2 density derived temperatures, and ail the 
temperatures have increased about 120°K from the undisturbed condi- 
tion in this region. Above that latitude the 6300 A temperature and the 
N_ density derived temperature increase at the rate of about 8°K per 
degree magnetic latitude. The total density derived temperature for 
higher latitudes has a slight increase at the rate of about 2°K per 
degree magnetic latitude. The increase in temperature over the 
undisturbed magnetic doncition at 70° magnetic latitude is about 350°K 
for the 6300 A temperatures and the N 2 density derived temperature, 
and only a few degrees for the total density derived temperature. The 
explanation for the total density derived temperature will be made 
evident when the variations of and O densities are discussed in 
Section 3,6. The 6300 A temperature and the density derived 
temperature disagree above 75° magnetic latitude where the 6300 A 
temperature continues to increase, whereas the density derived 
temperature decreases. 

Figure 3 shov/s the magnetic latitude variation of the 63 00 A 
temperature measurements and the N 2 and total density derived 
temperatures during the large disturbance caused by the great storm 
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of 8 March 1970 when the average value of Kp was 8o 15 and the 
maximum, value was 9.0, The 6300 A temperatures and the 
density derived temperature exhibit close agreement except from 
25° to 40° magnetic latitude where the density derived tempera- 
ture decreases slightly and then increases sharply at about 170°K per 
degree magnetic latitude while the 6300 A temperature increases at a 
relatively constant rate of about 26°K per degree magnetic latitude, 
and below 10° magnetic latitude where the difference is about the 
same as for the undisturbed condition. 

At this very high level of magnetic activity, it is apparent 
that the temperatures derived from total density fail to show the 700 K 
increase between the equator and the pole. The total density derived 
temperatures also show a large negative gradient in the 30° to 45° 
magnetic latitude region which is opposite to the variation shown by 
the 6300 A temperatures, 

3 , 2 Variations of Temperatures as a Function of Magnetic Index Kp 

Figure 4 shows the rate of change of temperature with Kp for 
the 6300 A temperatures and the temperatures derived from the 
and total densities. They were obtained using 11 days, June 14 to 24, 
1970. During this period the average value of Kp was 2.2 and the 
maximum value was 5.3. The AT/ A Kp values were obtained by 
least square fitting the data to the equation 

T = A + B SjQ ^ + (AT/AKp)Kp (4) 

where ^ is the 10,7 cm solar fl\ax index. 
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The coefficients of the least squares fits to the equation 

AT/AKp = K^+Kjsin0 (5) 

where 9 is the magnetic latitude, 
are given in Table 1 for comparison with the Roemer (1972) and 
Jacchia (1971) coefficients » It is apparent that the 6300 A temperature 
coefficient is much larger than that given by Roemer (1972) or 
Jacchia (1971), The total density temperature coefficient Kj is an 
order of magnitude smalle'r than that for the 6300 A and N 2 derived 
temperature coefficients and in better agreement with the Jacchia 
(1971) than the coefficients from Roemer (1972), 

Table 1: AT/AKp Coefficients for Summer 1970 Afternoon 



«0 


6300A 

-3,5 

38,6 

^2 

-6,2 

31,2 

Total Density 

7,3 

3,8 

Roemer (1972) 

17,9 

21,4 

Jacchia (1971) 

28, 0 

0,0 

Figure 5 shows 

the rate of change 

of temperature with Kp 


versus magnetic latitude in the same manner as Figure 4. These were 
obtained using data from 11 days, March 1 to 10, 1970, when the 
average value of Kp was 3,0 and the maximum value was 9o 0, The 
local time was around 800 hours. Table 2 gives the Kq and K^ 
coefficients for this period obtained in the same manner as Table 1, 
The AT/ AKp for the 6300 A and the derived temperatures are 
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Figure 5 AT/AKp versus magnetic latitude derived 
from 6300 A temperatures, and total 
density derived temperatures, Jacchia (1971) 
and Roemer (1972), for 1970 March AM 
local time 
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considerably larger than the total density derived temperature cases 
of Roemer (1972) and Jacchia (1971). The total density AT/ AKp is 
again closer to Jacchia (1971) than Roemer (1972). While the agree- 
ment between the 6300 A and N 2 density derived temperature is an 
order of magnitude better than for the total density derived tempera- 
ture AT/ AKp, 

Table 2; AT/ AKp Coefficients for March 1970 Morning 



Kq 

^1 

6300 A 

20.45 

54, 95 

N., 

-3.03 

101, 95 

2 

Total Density 

8.01 

11.18 

Roemer (1972) 

17.9 

21.4 

Jacchia (1971) 

28.0 

0, 0 


it appears that in this period the density derived temperature over- 
estimates the effect of magnetic activity at high magnetic latitude and 
underestimates it at low magnetic latitude. 

The March morning coefficients Kq and Kj are larger than the 
June afternoon coefficients. To examine the various factors influ- 
encing these coefficients, other periods were analyzed using the 
6300 A temperatures and the results are shown in Figure 6 and the 
coefficients listed in Table 3. A suggestion by Blamont and Luton 
(1972) that the effect of magnetic activity is stronger in the morning 
than in the afternoon is confirmed by both the K^ and K^ coefficients 
for the Fall 69 and the Summer 70 cases, A small seasonal variation 
can be detected in the Kq coefficient but it is not apparent in the K^^ 
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Figure 6 AT/AKp versus magnetic latitude calculated 
from 6300 A temperatures for Fall 69, 
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coefficients which would indicate that the global effect may be weaker 
in the summer but that the latitudinal variation is not changing with 
the seasons. Although the local time variation is clear, the seasonal 
results are not decisive, and much morns data should be investigated 
for other periods before hard conclusions can be drawn. 

Table 3: AT/ AKp Coefficients for Fall 69, Spring 70 and 
Summer 70 with AM and PM Local Times 



Kq 


Fall 69 AM 

13.6 

35.41 

Fall 69 PM 

11.76 

27. 88 

Spring 70 AM 

20,45 

54. 95 

Summer 70 AM 

o 

00 

59. 13 

Summer 70 PM 

- 3.5 

38. 57 


3.3 Fourier An alysis of Quiet Day Latitude Variation at Solstice 
and Equinox 

The OGO-6 satellite 6300 A temperature measurements were 
made during 3 solstices and 2 equinoxes. The orbit v/as such that a 
global distribution of morning and evening temperatures were obtained 
on each occasion. Thus it is possible to observe a global latitude 
temperature variations for each solstice and equinox, ana to determine 
the coefficients of a third order fourier series which best fit the 
global temperatures. 

The temperature data were sorted and averaged in 20 latitude 
ranges v/hich overlapped each other by 10 latitude, and a third order 
fourier series centered on the north pole was fitted to the averages. 
The results of the analysis are presented along with the averaged data 
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points in Figures 7, 8, 9, 10 and 11, and the fourier coefficients in 
Table 4o Due to a smaller number of data points, the Summer 69 
and Winter 69 periods required several days to obtain a global distri- 
bution of temperatures 0 

The fourier series analysis will allow others to compare 6300 A 
undisturbed global temperature variations with model temperatures 
and other experimentally measured temperatures, 

3 , 4 Anomalous. Temperature Behavior in the Southern Hemisjgher e_ 

Figures IZ, 13, 14 and 15 show the longitude variation of the 
6300 A morning temperature for 10° magnetic latitude ranges in the 
southern hemisphere for the interval 14-22 March 70 when the average 
value of Kp was 1,07 and the maximum value of Kp was 2„ 7 , The 
figures indicate an area from -20 to -40 magnetic latitude and 10 to 
60° west longitude which is consistently 300 K to 900 K hotter than 
the rest of the data for the period. This area of enhanced tempera- 
tures corresponds to the South Atlantic magnetic field anoinallyo 
Further study into other periods of low magnetic activity showed that 
the same area exhibited significantly higher temperatures only in the 
morning hours, usually before 1000 hours local timeo 

The densities and O densities were then investigated to find 
out if there was an N„ enhancement and O depletion corresponding to 
the same magnetic latitude and longitude o The results showed no 
significant variation of either the O densities at that location. 

It would seem that such large temperature gradients could not be 
maintained for long without producing a wind system which would 
largely equalize the pressure, ft must be tenatively concluded that 
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Table 4; Coefficients of Fourier Series Fitted to Quiet Day Orbital 
Temperature Variations at Equinoxes and Solstices 


T = 

Aq + sin 0 

+ Aj cos 0 

+ cos 

2 0 + sin 

26 + A 3 

cos 3 0 + 

cos 3 0 


■^0 



^2 

»2 

^3 


SjQ ^ 127.4 
Summer 69 

1019.48 

103. 92 

115, 79 

12. 56 

31. 82 

-0. 93 

-37.93 

Sio= 144.5 
Fall 69 

1036, 86 

-17. 54 

108, 00 

71.84 

-14. 87 

- 19 . 05 

-38. 93 

SjQ = 158.7 
Winter 69 

992, 89 

-105. 77 

95.55 

-61. 59 

-87, 09 

-18,91 
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S,o= 142.2 
Spring 70 

1005, 39 

2.41 
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JLU 

Summer 70 

o^ 

00 
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63 00 A temperature data versus longitude for magnetic latitude 
between -50° and -40° 
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the observed teraperatures do not correspond to the kinetic tempera- 
ture of the ambient gaSo Before the source of the emissions can be 
traced to its origin, a detailed examination of the spectra, should be 
made. This project could be of great interest and further study 
appears to be warranted, 

3 , 5 Temperatures at High Geomagnetic Latitudes 

Biamont and Luton (197Z) hav<» discussed the evolution of the 
temperatures at high magnetic latitudes during the period from 
September 23, 1969 to October 7, 1969 which included two severe 
geomagnetic storms. In this section an attempt will be made to 
obtain some average characteristics of the polar thermosphere under 
different seasonal and magnetic conditions. 

At high geomagnetic latitudes it has been shown that 
geophysical effects are represented most directly in a coordinate 
system based on geomagnetic latitude and geomagnetic time. In the 
present work the method of Agy (1965) v/as used as described in 
Appendix I of Matsushita and Campbell (1967) to calculate these param- 
eters. 

Temperatures at high magnetic latitudes are very variable and 
care must be taken in using satellite data to obtain time variations. 

The satellite precesses through 8 minutes in local time each day and 
storm conditions can cause an enhancement in temperature at the 
particular geomagnetic time when these measurements were made. 

To minimize, as far as possible, these fluctuations, it was decided to 
average the data over both hemispheres and quite extensive time 
periods , 
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The data way subdivided into ranges of geomagnetic latitude at 
5° intervals from 60° to 90° and in intervals of one hour in geomag- 
netic time, and in Kp ranges of 0 to 1» 3, 1„3 to 4 and 4 to 6c9« The 
small number of data points above 85° caused these temperatures to 
fluctuate greatly and thus measurements from 85° to 90° were 
averaged in 6 hour intervals „ 

The summer solstice data are shown in Figures 16, 17 and 18<. 
The days used for the northern hemisphere were 6 June 1969 to 
5 August 1969 and 7 May 1970 to 5 August 1970, and the days used to 
get the southern hemisphere summer were 5 November 1969 to 
3 February 1970o The tabulated values are given in Tables 5, 6 and 
7„ The accuracies shown are the estimated accuracy of each measure- 
ment divided by the square root of the number of measurements. 

The equinox data are shown in Figures 19> 20 and 21, The 
days used for these figures were 6 August 1969 to 4 November 1969, 
and days 4 February 1970 to 5 May 1970, The tabulated values are 
given in Tables 8, 9 and 10, and the accuracies are computed in the 
same way as the summer solstice. 

The winter solstice data are shown in Figures 22, 23 and 24, 
The days used are the same as the summer solstice case with the 
appropriate polar region selected. The tabulated values are given in 
Tables 11, 12 and 13 with the accuracies computed in the same way as 
the summer solstice. 

The seasonal effect at high geomagnetic latitudes appears to be 
greatest for low Kp values. The noon sector temperatures for 0 to 1,3 
Kp during summer were about 100°K hotter than winter, and during 
equinox were about 75°K hotter than winter. The midnight sector 
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Figure 16 Exospheric temperature as a function of 
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Exospheric temperature as a function 
of geomagnetic latitude and geomagnetic 
time for Kp range of 1.3 to 4.0 in Summer 
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Figure 18 Exospheric temperature as a function of 

geomagnetic latitude and geomagnetic time 
for Kp range of 4. 0 to 6. 7 in Summer 
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Table 6: Exospheric Temperatures and Accuracy of Averaged Data Tabulated by 
Geomagnetic Latitude and Geomagnetic Time for Summer with Kp Rang 
of lo3 to 4„ 0 


Geomagnetic Time Hours 

St^deT'" 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 


85 - 90 


80 - 85 


75 - 80 


70 - 75 


65 - 70 


60 - 65 


1174 1189 1232 1171 1278 1315 1340 1202 1072 1369 1085 1126 1145 1319 1220 1256 1166 1212 1233 1318 1335 1299 1^4 1214 

485 ±42 ±604120 ±69 ±49 ±49 ±85 ±60 ±60 ±69 ±54 ±69 *60 ±69 *60 *60 *85 


-IS ‘1^1 ‘ill ‘S ‘SI ‘S”o ‘SI ‘S! ‘S? ‘SI ‘SI ‘Si ‘Si ‘SI ‘SI ‘SI ‘31 ‘31 ‘SI ‘Si ‘SI 


1374 0 1228 1142 1196 1178 1225 1213 1198 1198 1182 1247 1247 1189 1191 1214 1264 1271 1267 1241 1317 1306 m6 m2 

±69 ±25 ±28 ±24 *24 ±29 ±24 ±20 ±21 ±25 ±30 ±31 ±26 ±26 ±26 ±38 ±31 ±36 ±40 ±30 ±3b ±45 ±45 


0 1302 1224 1243 1205 1214 1226 1227 1150 1141 1118 1109 1134 1160 1226 1242 1239 1208 1290 1294 12J4 1277 mi 
«4 ±21 ±24 ±21 ±20 ±17 ±22 ±21 ±21 ±30 ±28 ±26 ±24 ±23 ±22 ±24 ±25 ±24 ±26 ±26 460 ±85 


0 1222 1272 1182 1203 1177 1155 1141 1115 1066 1064 1165 1167 1247 1271 1260 1251 1247 1259 1^32 1^9 
0 As ±ll il8 fll ±14 ±23 ±24 ±20 ±28 ±29 *23, ±29 *20 ±23 *17 ±19 *21 *24 ±30 ±60 


0 1322 1269 1239 1211 1187 1175 1113 1101 1089 1071 1159 1189 H64 1222 1227 1217 1238 1273 im 
0 ±85 ±15 ±12 ±12 ±16 ±22 ±25 ±24 ±22 ±24 ±22 ±24 ±24 ±18 ±16 ±14 ±17 ±19 *69 
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Geomagnetic 
Latitudes 1 


Geomagnetic Time Hours 
10 11 12 13 14 


19 20 


22 23 


0 1178 1326 
0 ±120 ±120 


0 1364 
0 ±120 


0 1072 
0 ±120 


0 1427 
0 ±120 


0 1494 
0 ±120 


1552 1648 1571 1336 1640 1325 1340 1250 1291 1332 1330 0 1459 

±120 ±120 ±85 ±69 ±85 ±60 ±60 ±85 ±120 ±69 ±69 0 ±120 


0 1340 1262 1258 
0 ±120 ±120 ±85 


0 0 0 1163 0 0 

0 0 0 ±120 0 0 


0 1173 1440 1394 1248 1295 1353 1414 1373 1456 1379 1362 1368 
0 ±85 ±85 ±42 ±85 ±54 ±60 ±54 ±38 ±69 ±42 ±54 ±49 


0 1186 1243 1292 
0 ±120 ±85 ±69 


0 1256 1368 
0 ±120 ±120 


0 0 1427 1297 1401 1337 1375 1322 1253 1446 1213 1315 

0 0 ±54 ±60 ±54 ±49 ±38 ±49 ±45 ±45 ±85 ±54 


0 1345 1203 1219 1346 
0 ±38 ±85 ±60 ±45 


0 1385 1315 
0 ±120 ±69 


0 1254 1376 1299 1331 1342 1360 1259 1216 1476 1137 1383 1348 1341 1308 1344 1319 1340 1306 1273 
0 ±120 ±49 ±49 ±30 ±40 ±85 ±45 ±60 ±69 ±60 ±69 ±54 ±54 ±69 ±35 ±85 ±69 ±60 ±i85 


0 1417 1340 1289 1259 1255 1195 1301 1325 1166 1274 1304 1356 1297 1354 1288 1305 1281 1308 
0 ±38 ±35 ±30 ±42 ±69 ±120 ±54 ±69 ±69 ±60 ±69 ±60 ±60 ±45 ±85 ±49 ±54 ±120 
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Figure 19 


Exospheric temperature as a function of 
geomagnetic latitude and geomagnetic time 
for Kp range of 0 to 1.3 at equinox 
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Exospheric temperature as a function of 
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for Kp range of 1.3 to 4. 0 at equinox 
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Figure 21 


Exospheric temperature as a function of 
geomagnetic latitude and geomagnetic time 
for Kp range of 4.0 to 6.7 at equinox 






Table 8; Exospheric Temperatures and Accuracy of Averaged Data Tabulated by 
Geomagnetic Latitude and Geomagnetic Time for Equinox with Kp Range 
of 0 to 1,3 


Geomagnetic Time Hours 

^2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 


85 - 90 


'ai 'a? ‘a? 'a? ‘a? ‘a? ‘ai aa ‘u? at ‘ai ‘ai -aj 'a? 


80-85 


10,3 10« 10,3 103. 1M3 .0^3 1110 110, m3 1083 1041 1030 1100 10,1 1136 1.36 im 1H8 1068 m3 


75 - 80 


‘a! 'a? *a^ ‘a? ‘a? ‘ai ’a *a^ ‘si *ai aa aa aa aa 'itt 'a? -a? *ia ‘a? ‘a? as as ‘at 
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0 1016 1022 1004 1012 978 1032 1010 1004 1033 988 986 1057 1088 1099 1081 1123 1072 1117 U30 1124 1«9 
0 ±49 ±33 ±24 ±18 ±15 ±10 ±10 ±10 ±10 ±11 ±15 ±14 ±17 ±13 ±15 ±15 ±16 ±16 ±17 ±ZO ±45 


60 - 65 


0 1083 996 982 965 968 1027 985 983 987 995 1023 1039 108^112 1059 1061 1121 1086 10^^ 896 

0 ±69 S 3 S 7 ±21 ±15 ±9 ±9 ±8 ±10 ±13 ±15 ±14 ±16 ±12 ±12 ±13 ±15 ±15 ±16 ±33 ±85 
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Latitudes 


id Accuracy of Averaged Data Tabulated by 
lomagnetic Time for Equinox with Kp Range 


s-amagnetic Time Hours 

0 11 12 13 14 15 16 17 18 19 20 21 22 23 24 


,82 1158 1147 1247 1248 1113 1232 1069 
t 36 ±29 ±36 ±40 ±33 ±49 ±30 ±40 


L 91 1175 1198 1190 1164 1125 1097 1141 
fc 21 ±22 ±18 ±21 ±18 ±17 ±19 ±18 


128 1159 1179 1170 1131 1091 1097 1160 
tl 3 ±12 ±12 ±15 ±13 ±12 ±14 ±16 


1123 
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±9 ±11 ±11 ±11 ±12 ±11 ±10 ±12 ±14 ±15 ±16 ±14 ±15 ±22 i 38 


048 1076 1047 1035 1097 1118 1160 1147 1102 1135 1148 1168 1107 1222 0 

±8 ±9 ±10 ±14 ±12 ±12 ±11 ±11 ±12 ±13 ±14 ±12 ±15 ■ 0 


061 1012 999 1022 1014 1046 1128 1103 1113 1139 1112 1078 1088 9 W 0 

±7 ±9 ±11 ±14 ±12 ±14 ±11 ±11 ±11 ±12 ±12 ±14 ±26 ±60 0 

















Table 10: Exospheric Temperatures and Accuracy of Averaged Data Tabulated by 
Geomagnetic Latitude and Geomagnetic Time for Equinox with Kp Range 
of 4.0 to 6.7 
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Geomagnetic Time Hours 
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16 
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20 
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22 

23 
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0 
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0 1263 
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Figure 22 Exospheric temperature as a function of 

geomagnetic latitude and geomagnetic time 
for Kp range of 0 to 1 . 3 in Winter 
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Figure Z3 Exospheric t-emperature as a function of 

geomagnetic latitude and geomagnetic time 
for Kp range of 1.3 to 4.0 in Winter 
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Exospheric temperature as a function of 
geomagnetic latitude and geomagnetic time 
for Kp range of 4. 0 to 6. 7 in Winter 







Table 11: Exospheric Temp era tiires and Accuracy of Averaged Data Ta,bulated by 
Geomagnetic Latitude and Geomagnetic Time for Winter with Kp Range 
of 0 to 1 . 3 


Geomagnetic 
Latitudes 1 

2 

3 

4 

5 

6 

7 

8 

Geomagnetic Time Hours 
9 10 11 12 13 14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

85 - 

90 

795 

0 

970 

0 

830 

1023 

865 

1282 

792 

854 

936 

928 

798 

1034 

901 

921 

1030 

972 

741 

1114 

0 

937 

0 

0 



±120 

0 

±120 

0 

±120 

±69 ±120 

±60 ±120 

±85 ±120 ±120 

±85 

±69 

±85 

±85 

±54 

±69 ±120 

±60 

0 

±85 

0 

0 

80 - 

85 

1277 

897 
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730 
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Table 12: Exospheric Temperatures and Accuracy of Averaged Data Tabulated by 
Geomagnetic Latitude and Geomagnetic Time for Winter with Kp Range 
of 1.3 to 4. 0 


Geomagnetic 

Latitudes 


85 - 90 


80 - 85 


75 - 80 


70 - 75 


65 - 70 


60 - 65 


1 

2 

3 

4 

5 

6 

7 

8 

Geomagnetic Time Hours 
9 10 11 12 13 14 15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

1532 1369 

0 

0 

879 

0 

1254 1149 

1040 

947 983 1257 

1068 1008 1598 1157 

1100 

1036 

876 1088 

882 

850 

1056 

1116 

±120 

±85 

0 

0 

±60 

0 

±69 

±60 

±49 

±85 ±60 ±120 

±54 

±85 ±120 

±49 

±49 

±60 

±85 

±42 

±85 ±120 ±120 ±120 

1085 

0 

1006 

924 

936 1170 

1052 1026 

1052 1046 1116 1043 

1070 

1127 1073 

1043 

1034 1037 

1051 

1117 

1081 

1006 

0 

732 

±85 

0 

±85 

±85 

±49 

±54 

±42 

±35 

±31 

±40 ±26 ±23 

±30 

±26 ±33 

±28 

±38 

±33 

±49 

±60 

±85 

±69 

0 ±120 

857 

0 

851 

889 

1044 

995 

1072 

1068 

1087 

1093 1093 1093 

1025 1081 1096 1113 

1037 

1023 

1056 1036 

1010 

693 

0 

0 

±85 

0 

±60 

±54 

±42 

±29 

±30 

±25 

±22 

±20 ±22 ±20 

±17 

±20 ±21 

±22 

±20 

: t 26 

±42 

±54 ±120 ±120 

0 

0 

0 

0 

0 

942 

921 1010 

961 

1008 

1036 

1029 1036 1041 

1018 

1070 1090 

1091 

1034 

1010 

1001 

956 

0 

0 

0 

0 

0 

0 

0 

±85 

±49 

±36 

±33 

±18 

±22 

±17 ±18 ±16 

±15 

±15 ±16 

±17 

±14 

±16 

±31 

±69 

0 

0 

0 

0 

0 

0 

853 

800 

934 1091 

992 1016 

993 

984 1011 938 

988 1042 1053 

1048 1027 

1020 

1011 

879 

807 

0 

0 

0 

0 

0 

±120 

±85 

±45 

±69 

±32 

±19 

±18 

±18 ±16 ±14 

±13 

±14 ±15 

±15 

±15 

±12 

±35 

±45 ±120 

0 

0 

0 

0 

0 

0 

812 

882 

1131 

934 

981 

950 1027 948 977 

935 

1007 965 1007 

1016 

1022 

941 

833 

0 

0 

0 

0 

0 

0 

0 

±69 ±120 

±85 

±42 

±22 

±24 

±18 ±15 ±13 

±10 

±12 ±12 

±15 

±16 

±14 

±40 

±85 

0 

0 

0 

0 





Table 13: 


Exospheric Temperatures and Accuracy of Averaged Data Tabulated by 
Geomagnetic Latitude and Geomagnetic Time for Winter with Kp Range 
of 4, 0 to 6o 7 


Geomagnetic Time Hours 


22 23 24 



85 - 90 


1719 0 

il20 0 


0 0 1559 

0 0 ±120 


80 - 

85 

0 

0 

1530 

1202 

0 



0 

0 

±120 

±85 

0 

75 - 

■ 80 

0 

0 

863 

1192 

942 

0 

0 

±120 

±54 

±85 

70 - 

• 75 

0 

0 

1053 

1190 

1193 

0 

0 

±85 

±85 

±85 





105 



65 ■ 

- 70 

0 

0 

1058 

833 

916 



0 

0 

±120 

±120 

±85 

60 

- 65 

0 

0 

0 

867 

1448 

0 

0 

0 

±120 

±85 


796 1156 1299 1185 ViH libO luit uiu ---- 

±85 ±69 ±120 ±54’ ±69 ±69 ±85 ±60 ±69 ±69 ±60 ±120 ±120 ±85 ±120 


1270 1220 107Z 115± LlO! i-io-i xioa xjvu ---- ---- ---- 
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±60 ±49 ±54 ±60 


±60 ±45 ±38 ±42 ±45 ±38 ±40 ±40 ±40 


0 1270 
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0 0 
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temperatures for Kp from 0 to 1.3 during summer were about 325 
hotter than v/inter, and during equinox were about 225 K hotter than 
winter. The summer and equinox seasonal variation for the 5 to 6, 7 
range of Kp were about one- third as large as for the 0 to 1.3 Kp 
range. 

The effect of magnetic activity on the temperatures during 
summer showed an increase of 80 K for the 1,3 to 4 Kp range from 
the 0 to 1,3 Kp range temperatures, and an increase of 200°K for the 
4 to 6.7 Kp range from the 0 to 1.3 Kp range temperatures. The 
equinox temperatures showed increases of 75°K for the 1,3 to 4 Kp 
range and 250°K for the 4 to 6, 7 Kp range from the 0 to 1 . 3 Kp range 
temperatures. The winter temperatures showed increases of 100°K 
for the 1,3 to 4 Kp range and 350°K for the 4 to 6.7 Kp range from 
the 0 to 1.3 Kp range temperatures. 

The geomagnetic local time effect is evident from the large 
maximum of temperature around 1900 to 2100 hours and a small maxi- 
mum around 700 to 1000 hours. The noon sector is generally the 
minimum temperature area. The absolute maximum temperature 
area can not be denoted as there is little data available in the geomag- 
netic local times near 100 hours, however, there are indica,tions that 
the maximum is in the pre-midnight sector. Caution is necessary 
regarding the geomagnetic local time effects as the separation of 
temporal and magnetic activity effects is difficult. The 2100 hour 
maximum and the noon minimum temperatures appear consistently 
throughout the data. The 900 hour maximum temperature appears 
often, but as the variation is smaller, it is not as apparent as the 


2100 hour ma.ximum. 
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Figure 25 shows the incidence of auroral forms in 1964 and 
1965 and modified to give a resolution similar to the temperature data„ 
This can be used to compare with the polar maps already presentedo 
The outstanding feature of Figure 25 is that the incidence of auroral 
forms is highest near 2100 hours and 800 hours geomagnetic local 
time which correspond to the aforementioned geomagnetic local time 
temperature maxima„ The behavior in general is, however, very 
different^ The exospheric temperatures do not decrease within the 
auroral oval and show much less spatial structureo In almost all the 
graphs it can be seen that the magnetic pole is at a temperature close 
to the maximum in the region above 60 „ It would appear that the 
effects of convection and conductivity are important enough to 
redistribute the energy deposited throughout the polar region and 
particularly, within the auroral ovaL The gradients are quite large, 
however, compared with those found elsewhere,, Even under quiet 
conditions with Kp lo3 they are as large as 200 K in 1500 km in the 


noon sec tor o 
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0 hrs 


Figure 25 Incidence of auroral forms during 1964 and 
1965 after Stringer and Belon (1966) 
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3,6 Description of the Temperature, Density, O Density, 

Pressure and Boundary Oxygen Density Variation During the 

Great Storm of 8 March 1970 

Figures 1, 2 and 3 in Section 3, 1 have already given an indica- 
tion of the variation of the temperature in the northern hemisphere for 
the period around 8 March 1970, Figure 26 shov/s the 63 CA tempera- 
ture from -60° magnetic latitude to 90° magnetic latitude for the 
undisturbed day 11 March 1970, the moderately disturbed days 6 and 

7 March 1970, and during the peak temperatures on 8 March 1970, 

The large temperature increase of 600°K to 700°K is evident at 
both north and south high latitude regions. The increase at the 
equator is only about 200°K, The latitudinal temperature gradient for 
the moderately disturbed day is about 7°K per degree latitude to the 
pole and slightly larger than the undisturbed day. The gradient for 
greatly disturbed day is about 11°K per degree latitude from 10° to 
50° magnetic latitude and then fluctuates near 0°K per degree latitude 
from 50° to 90° magnetic latitude. This indicates that the energy is 
being deposited above 50° magnetic latitude and then transported to 
the lower latitudes. 

Figures 27, 28 and 29 show a contoured view of the 6300 A 
temperature, density and O density versus magnetic latitude and 
Universal Time for five days during and after the great storm of 

8 March 1970, Day No, 67 is the greaTly disturbed day and day No, 71 
is an undisturbed day which may be used for comparison. Unfortu- 
nately, the density data was lost for day No, 67 until U,T„ of 1900 
hours. The Kp and AE indices for the interval are shown in Figure 30, 

The 63 00 A temperatures in Figure 27 illustrate the rapid 
response of the thermospliere to a large heat source. The time of 


X - II MARCH 70. 4 - 24hrs U.T. 
0-6-7 MARCH 70, 0-24hrs U.T. 
+ - 8 MARCH 70. 16-24 hrs U.T. 




U.T. 4 12 20 4 12 20 4 12 20 4 12 20 4 12 20 

!970 DAY# 67 68 69 70 71 

Figure 27 63 00 A tempera ture (in hundreds of °K) contour map 

representing the 300 km variations with magnetic 
latitude and universal time for five days in March 1970 
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Figure 29 Log,„ (0 number density (m )) contour map representing 
the 40Q Ian variations with magnetic latitude and universal 
time for five days in March 1970 
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response at the high latitudes is less than one hour while the response 
time at the equator is around three or four hours« These reduced 
time responses from previously published results (Roemer, 1972 and 
Alcayde, et al, , 1974) are certainly a result of the expansion of the 
auroral oval and the equatorward movement of the heat source, as was 
mentioned previously and can be seen in Figure 27, 

The 400 km density contours of Figure 28 also show au 
order of magnitude difference between the densities at the pole and the 
equator. The latitude variation on day No. 71 is very much smaller. 
The general variation of the density contours can be seen to be 
quite similar to the 6300 A temperature contours, though there are 
disagreements at several locations as had been previously pointed out 
in Sections 3,1 and 3.2 when the N2 density derived temperatures 
were compared to the 6300 A temperatures. 

The atomic oxygen 400 km density contours are presented in 
Figure 29 and can be seen to have variations v/hich are strikingly 
different from the density and 6300 A temperature contours. The 
most important feature of Figure 29 is the moderate increase of 
atomic oxygen below 35° magnetic latitude and the large decrease of 
atomic oxygen above 3 5° magnetic latitude. Thus, the variation of 
atomic oxygen with moderate or large magnetic activity can be 
increasing or decreasing depending on the magnetic latitude which is 
being investigated. This behavior of the atomic oxygen explains the 
variations of the total density derived temperatures discussed in 
Sections 3, 1 and 3,2 and also affects the pressure variations. 

The pressure during undisturbed, moderately disturbed and 
greatly disturbed March 70 morning periods is shown in Figure 31, 


PRESSURE (ntm 


X - II MARCH 70 
4-24hrs U.T. 

□ - 6-7 MARCH 70 
0-24hrs U.T 

+ - 8 MARCH 70 
16-24 hrs U.T 


10 ® L 
- 10 ® 


20® 30® 40® 50® 60® TCP 

MAGNETIC LATITUDE 

7.7 8.0 8,5 

LOCAL TIME 


Figure 31 Variation of pressure at 400 km versus 
magnetic Latitude for 'mdis terbed, 
moderately disturbed and gr ea.tly dis turbed 
period in. March 1970 
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The undisturbed and moderately disturbed periods indicate a low 
pressure region around 20° to 30° magnetic latitude and higher 
pressures towards the poles „ The greatly disturbed pressure varia- 
tion is generally flatter with large gradients at several locations. 

The pressure appears to equalize very quickly during large disturb- 
ances which would indicate large but short-lived winds are generated 
and the atomic oxygen quickly transported to low pressure areas to 
equalize the global pressure. 

Figure 32 shows the undisturbed^ moderately disturbed and 
greatly disturbed magnetic latitude variations of the atomic oxygen 
near 120 km. The boundary atomic oxygen density was calculated 
using the following equation 

n Z 

n(OQ) = n(0) (T/Tq) exp ( J m^g/kT dZ) (6) 

v</ 3re r e 

si(Oq' and T^ are the boundary density and temperature 
n.(0) and T are the 400 km density and temperature 
and the exponential term is evaluated using the N.^ density and 6300 A 
temperatures. Figure 32 shows a moderate increase in the boundary 
oxygen density at magnetic latitudes below 40° and a great depletion 
above 40° magnetic latitude. This would indicate winds up and away 
from the high magnetic latitudes down and into the low magnetic 
latitudes during magnetic disturbances. 

Figure 33 shov/s the 400 km atomic oxygen variations corres- 
ponding to the same intervals as Figure 22, It must be noted that the 
depletion at the high latitudes is smaller at 400 km and that the 
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Figure 32 Variation of lower boundary atomic 
oxygen versus magnetic latitude for 
undisturbed, moderately disturbed, 
and greatly disturbed periods in 
March 1970 
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Figure 33 Atomic oxygen at 400 km versus 
magnetic latitude for undisturbed, 
moderately disturbed and greatly 
disturbed periods, and 3 hours after 
the greatly disturbed period in 
March 1970 
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enhancement at the low latitudes is larger than the atomic oxygen 
near iZO km„ This is consistent with a transport of atomic oxygen up 
and away from the high latitudes and a subsidence at the low latitudes = 

3 o 7 Time Constants of Thermospheric Effects Due to Magnetic 
Activity 

In a magnetic storm, energy deposition occurs in the high 
latitude regions. While the AE and Kp indices provide a measure of 
the geomagnetic activity, it must be presumed that the effects will 
demonstrate both time delay and energy storage phenomena. Time 
delays result from the time necessary to transport the energy to 
lower latitudes by winds and gravity waves, and the storage effects 
will be associated with the time required to change the temperature 
gradients at lov/er altitudes. 

If it is assumed that the thermosphere has an integrating 
effect with regard to the magnetic activity, then the thermosphere 
would be responsive to the magnetic activity over some proceeding 
time interval, and. after an impulse of magnetic activity heating would 
relax to its original condition with some relaxation time. To test this 
hypothesis and attempt to determine the time constant which would 
relate the thermospheric response to magnetic activity, three 
different periods of activity were selected. The three intervals were; 
(1) Fall 69, 1 to 11 September 1969, (2) Spring 70, 1 to 11 March 
1970 and (3) Summer 70, 14 to Z4 June 1970, 

It has been shown in Section 3,4 that the initial response to a 
magnetic disturbance can be very rapid and that the increase in 
temperature can be related to the magnetic indices using delays 
ranging from less than an hour at 60° magnetic latitude to about 
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10 hours at the equator » If the period after the storm is observed, it 
is apparent that the temperatures do not return to the quiet values as 
quickly as the magnetic indices because of the energy storage in the 
thermosphere„ This is the cause of two common errors in dealing 
with the thermospheric temperatureso First, the model tempera- 
tures based on a one-to-one relation to Kp will underestimate the 
temperatures following a storm. Second, experimenters may use 
quiet day data to construct models of the thermosphere from periods 
v/here the Kp was low, but with temperatures that were enhanced 
due to the relaxation time constant effect of a preceeding storm 
period. 


It is apparent from the data that what is needed is a method of 
using the magnetic indices so that the beginning of the sto^m is well 
defined and the period after influenced by the relaxation time effect. 
With this in mind, a new index has been employed to analyze the 
6300 A temperatures, A factor W was introduced to account for the 
magnetic variation where W is evaluated by 

n n 


W(t, 


^d’ ~ y (-31 / t) / ^ exp (-31 /t) 


i=0 


i=0 


(7) 


where 


t = time (U,To), t^ = time delay (hrs), t = relaxation time 

constant (hrs), and Kp(t-t^-3i) is the appropriate magnetic index 

n 

weighted by exp (-31/t) and ^ exp (-Si/r) is a normalizing factor so 

i=0 

that W values for different t's are of comparable magnitude. 
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Figure 34 shows the error of the fit for the different periods 
already specified as a function of the relaxation time constant t« 

T = 0 is just the variation with Kp, A linear function of W was used in 
this study. The best fit for Spring and Fall corresponds to t = 6 hrs, 
while the best fit for Summer was t = 12 hrs. There appears to be a 
definite seasonal variation in the relaxation time constant although 
there doesn’t seem to be a local time variation. 

Figure 35 shows the magnetic latitude variation of the 
at/ aw coefficient for Summer 70 AM and PM, As expected, the 
coefficients for 700 hrs local time are significantly larger than those 
for 1800 hrs and the high latitude coefficient larger than the Low 
latitudes , 

Figure 36 shows the magnetic latiti’de variation of the 6300 A 
temperature with W = 0 for the Summer 70 period using the fitting 
function, A comparison of Figure 36 with Figure 11 shows an 
excellent agreement between the quiet day 23 June 70 temperature 
data and the results of the fit with W = 0, This indicates a good 
representation of the magnetic activity by the W factor, 

Som.e important conclusions can be drawn about the relaxation 
time of the thermosphere for magnetic activity. It does not appear 
that these times are longer than 12 hours and there does not seem to 
be a large latitudinal gradient in the relaxation time. 


STANDARD DEV. - ERROR OF FIT («K) 
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Standard deviation of least squares fit of 
6300 A temperature data versus relaxation 
time constant T for Fall 69, Spring 70 and 
Summer 70 for AM and PM local times 
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Figure 36 Exospheric temperature versus magnetic latitude for W = 0 and 
lOo 7 cm solar flux = 150 for Summer AM and PM local times 
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CHAPTER IV 
CONCLUSIONS 

4 , 1 Cotnparison of the 6300A Temperature with. Tem.peratu.res 
Derived from Total Density and from N^ Density 

It has been shown that the total density when used with static 
diffusion models with constant boundary conditions and shape param- 
eters gives an exceedingly poor mea.surement of the global variation 
of the neutral temperature under both quiet and disturbed conditions. 
Under quiet conditions the temperatures derived from total density 
show a variation from equator to pole of about 80 K compared with 
the equator pole variation of 300 K for the 6300 A temperatures and 
150°K for the temperatures derived from N^ density. Under disturbed 
conditions the temperatures derived from total density fail to describe 
the global variation of temperature. The depletion of atomic oxygen 
in the high latitudes ca.uses the temperatures derived from the total 
densities to be quite different from the kinetic temperatures as shown 
in Figure 3. 

The agreement between the 6300 A temperatures and the 
temperatures derived from N 21 densities was much closer than the 
total density derived tempejrature, but there were still major 
discrepancies. In agreement with incoherent scatter measurements 
at mid- latitudes, the N^, derived temperature was found to be 100°K 
to 200°K too high under quiet conditions. Daring disturbed conditions, 
the temperatures derived from the N^ densities agreed quite well with 
the 63 00 A temperatures, but the latitudinal temperature gradient in 
the Np^ derived temperature gradient was much greater in the region 
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below the auroral zone and the derived temperature poleward of 
the auroral zone decreased rather than increased as the 6300 A 
temperatures indicated. The reason for these discrepancies are 
probably more closely associated with the temperature gradient in the 
lower thermosphere rather than with transport effects. It is known 
that large changes in the temperature profile and total density occur 
in the region of IZO km due to the semi-diurnal tide. Large changes 
also appear to occur during conditions of high magnetic activity so 
that it is not surprising that there is not a one-to-one correspondence 
between exospheric temperature and 400 km density. The temperature 
derived from the N„ densities would be expected to give a closer 

C-t 

measurement of the temperature than any other gas because it is the 
major constituent in the region where eddy diffusion and transport are 
important, and because its mass density is close to that oi the mean 
molecular mass in this region, 

4 , 2 Temperature Variations as a Function of Kp Magnetic Index 
The results of the analysis of 6300 A temperatures and 
temperatures derived from density and total mass density are 
summarized in Tables 1, 2 and 3, 

The results shown in Table 3 clearly indicate a strong local 
time variation of the global Kq and latitudinal coefficients of 
AT/ AKp for both Fall 69 and Summer 70 storm periods. There is a 
moderate seasonal variation in the global Kq coefficient but no clear 
seasonal variation in the latitudinal coefficient. 

Tables 1 and 2 show that the latitude gradient in the 6300 A 
temperature enhancement is an order of magnitude larger than 
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indicated by the method used to derive the coefficients for Jacchia 
(1971)„ There was general agreement between the total density 
measurements and the total densities of Jacchia (1971), however, 
several important departures in detail were noticed. At the peak of 
the storm the total density decreased rapidly with increasing latitude 
in the region of the midlatitude trough. The reasons for these 
discrepancies were studied. It is apparent from Figure 32 that the 
changes are due to changes in the lower boundary conditions of atomic 
oxygen and represent the effects of transport. 

The 6300 A temperature AT/ AKp coefficients showed good 
agreement with the variations calculated using the temperatures 
derived from the N 2 density except that the density derived 
temperature variations underestimated the effects at low latitudes and 
overestimated the effects at high latitudes. 

4, 3 Equinox and Solstice Global Variations for Quiet Magnetic and 
Solar Conditions 

The Fourier series analysis of the quiet day global tempera- 
ture distribution allows an easy means of comparing the 6300 A 
temperatures measured during solstice and equinox conditions’ with 
model temperatures or other experimental temperature data. The 
coefficients are presented in Table 4, 

4, 4 Angmalous Southern Hemisphere Temperature Longitude 
Variation 

A large enhancement of the 6300 A temperature measurements 
in the area of the South Atlantic magnetic field anomally was demon- 
strated by Figures 12, 13, 14 and 15, The enhancement was of the 
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order of 1000°K„ No evidence of a corresponding variation of the 
O or N densities was seen<, Further investigation into the cause of 
this temperature effect is neededo 


A. ■; Temperature V ariation as a Function of. Ge omagneticJTHK^^ 
Gi somagnetic Latitude 

The results of the analysis of the 6300 A airglow temperatures 
as a function of geomagnetic time and geomagnetic latitude are shown 
in Section 3 A. The local time variation should be regarded with 
caution because of difficulties in separating the temporal and magnetic 
activity effects o The enhanced temperatures in the pre-midmght 
sector and the low temperatures in the noon sector seem to be real 
geophysical effects and are clearly seen in the summer and equinos: 
datao Latitudinal gradients are not large and high temperatures are 
maintained up to the magnetic polco Season variations are larger 
under low magnetic activity than for high magnetic activityo 

d 6 Tp.yrmeratu re„ Density. Pressure aEid_B ounda££,0^ 

and Afte_r_a Jja£geJvl^^ 

The great storm of 8 March 1970 was analyzed in detail in 
Section 3=6„ Figures 27, 28 and 29 showed the Universal Time and 
magnetic latitude variations of the temperature, density and O 
density at 400 km corresponding to the magnetic activity indicated in 


the magnetic indices shown in Figure 30o A fast response at all 
magnetic latitudes was shown to be on the order of one hour at 60 
magnetic latitude to about three hours at the equator « The cause of 
the short response time may be the expansion of the auroral oval to 
midlatitudes which would advance the source of heating closer to tne 


equator » 
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The 6300 A temperatures and N 2 densities demonstrated 
similar variations during the interval investigatedo The densities 
exhibited several sharp latitude gradients which were not present in 
the temperatures. 

The atomic oxygen variation was strikingly different from 
either the N 2 density or the 6300 A temperatures. Figures 32 and 33 
showed the boundary oxygen and 400 km oxygen density variations to 
be indicative of a transport mechanism which moves the atomic 
oxygen up and away from the high Ia.titudes towards the equator where- 
the movement subsides. 

The cause of the atomic oxygen behavior appears to be 
connected with the pressure. With the OGO-6 satellite,, independent 
measurements of the neutral temperature and density could be 
combined to obtain an excellent measure of the pressure. As seen in 
Figure 31 the north- south pressure gradient was generally small 
except for several large localized gradients at the height of the storm. 
It appears that the thermospheric pressure induces winds of such 
magnitu.de that the pressure gradients are minimized, 

4, 7 Temperature Relaxation Time Constant 

A new magnetic activity index W was introduced to allow the 
response of the thermosphere to magnetic storms to be modeled 
better. This index is defined in equation 7, It is based on the Kp 
index and it includes both a delay time and an exponential time delay. 

It was found that this index provided significantly better fits between 
temperatures and the index than did the Kp value itself. It was found 
that the relaxation time constant was of the order of 3 to 12 hours and 


not strongly latitude or local time dependent. The time delays varied 
from less than one hour at high latitudes to 10 hours at the magnetic 
equator , 
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